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Abstract 

The contribution of the low-lying nucleon resonances ^33(1232), Pii(1440) -Di3(1520) and 
5ii(1535) to the invariant mass spectra of di-electrons stemming from the exclusive processes 
pp — > pp e~^e~ and pn — > pn e~^e~ is investigated within a fully covariant and gauge invariant 
diagrammatical approach. We employ, within the one-boson exchange approximation, effective 
nucleon-meson interactions including the exchange mesons vr, rj, a, to and p as well as excitations 
and radiative decays of the above low-lying nucleon resonances. The total contribution of these 
resonances is dominant, however, bremsstrahlung processes in pp and, in particular, pn collisions 
at beam energies of 1 - 2 GeV are still significant in certain phase space regions. 
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I. INTRODUCTION 



The experimental study of di-electrons as penetrating probes in relativistic heavy-ion 
coUisions is aimed at identifying medium modifications of hadrons, in particular of the 
vector mesons p, cu and (f) [1]. Previous measurements of di-electrons in the reaction 
^^C -I- ^^C at kinetic beam energy of 1.04 AGeV performed by the DLS collaboration [2] 
have been confirmed recently by the HADES collaboration [3], at least in phase space 
regions covered by both experiments. Various transport models have been employed [4, 5] 
for understanding and interpreting the di-electron data [2, 3, 6]. Among the impor- 
tant sources for di-electrons in the low- mass region are tt'^, A and rj Dalitz decays and 
bremsstrahlung as well [4, 5] . 

The elementary cross section for virtual nucleon-nucleon bremsstrahlung with 7* — > 
e'^e~ as a subprocess in heavy-ion collisions was parameterized often within the soft- 
photon approximation (cf. [7]), which is appropriate at low kinetic energies, where the 
photon is quasi-real, but becomes questionably at higher energies and at higher virtual- 
ities of the 7*. Moreover, the soft- photon approximation preserves only approximately 
the gauge invariance, and the violation of gauge invariance increases with initial energy. 
In Ref. [8], based on previous investigations [9-13], a fully covariant and gauge invariant 
approach has been proposed to parameterize the bremsstrahlung amplitude in elementary 
pp and pn collisions. It was demonstrated that, in order to preserve the gauge invariance 
in pn reactions, one has to include additional diagrams with meson exchange currents and, 
for the couplings with field derivatives, to introduce contact terms, the so-called seagull 
or KroU-Rudermann [14] type diagrams. The resulting pn bremsstrahlung cross section 
was found to essentially differ from the one obtained within previous quasi-classical cal- 
culations. (This conclusion has been confirmed in Ref. [15].) The calculations reported 
in [4] utilized the bremsstrahlung cross sections of [8] and, indeed, are capable describ- 
ing perfectly the DLS [2] and the recent HADES di-electron data [3, 6] for the reaction 
i2q _|_ i2(-y Hence, one can assert that the so-called "DLS puzzle" originated from scarce 
knowledge of elementary cross sections used in transport models, in particular the ele- 
mentary nucleon-nucleon bremsstrahlung. 
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In covariant approaches, based on an effective meson-nucleon theory to calculate the 
bremsstrahlung of di-electrons from nucleon-nucleon scattering, the effective parameters 
have been adjusted to describe elastic nucleon-nucleon (NN) and inelastic NN NNn 
processes at intermediate energies. Excitations of resonances have been studied at the 
same time, and it is found that at intermediate energies the main contribution comes 
from A resonances (see also Ref. [16]), whereas excitations of higher mass resonances are 
often neglected. The role of higher mass and spin nucleon resonances at energies near the 
vector meson (p, cu and 0) production thresholds have been investigated for proton-proton 
collisions in several papers (see, e.g., Refs. [17, 18] and further references therein) with the 
conclusion that at threshold-near energies the inclusion of heavier resonances also leads 
to a good description of data. However, as demonstrated in Rcfs. [10, 17] calculations 
with a reasonable readjustment of the effective parameters can equally well describe the 
data without higher mass and spin resonances. In contrast, for di-electron production in 
photon and pion induced reactions, excitations of low-lying as well as heavier resonances 
can play a role [19]. 

In the present paper we investigate in some detail the role of nucleon resonances with 
masses close to the A for di-electron production in A'^A'^ collisions. Besides the A we 
consider the low-lying Pii(1440), Z^i3(1520) and 511(1535) resonances which are expected 
to contribute at larger values of the e+e" invariant mass and, therefore, can modify the 
shape of the e'^e" mass distribution at the kinematical limit. 

Our paper is organized as follows. In section II A we recall the kinematics and the 
general expressions for the cross section. The purely electromagnetic part of the cross 
section is considered in section II B, where the integration over the leptonic variable is 
performed analytically and an expression for the cross section is presented. In sections 
II C and II D, the effective Lagrangians and the problem of gauge invariance within the 
one-boson exchange model are discussed. Meson exchange diagrams and seagull terms 
are considered in this section as well. Results for the invariant-mass distribution of di- 
electrons stemming from pp and pn bremsstrahlung processes, where only nucleons and 
mesons are involved, arc reported in section II E. The role of resonances is investigated in 
section III. In section III A, the contribution of the A isobar is considered. In particular. 
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the choice of the couphng constants together with the off-mass shell parameters is dis- 
cussed. A comparison of the contributions from bremsstrahlung and A is presented also in 
this subsection. The nucleon resonances with spin 1/2 and 3/2 are considered in sections 
III B and III C, respectively. The adjustment of effective parameters to experimental data 
and the parametrization of the energy dependence of the resonance widths are reported 
in detail. The individual contributions of each resonance are analyzed. The total cross 
section as a coherent sum of bremsstrahlung and resonance contributions, including all 
interference effects, is presented for two experimentally relevant kinetic energies in pp and 
pn collision. The summary and conclusions can be found in section IV. 



II. DI-ELECTRONS FROM NN COLLISIONS 



A. Kinematics and Notation 

We consider the exclusive e^e^ production in NN reactions of the type 

Ni{Pi) + N2{P2) ^ N[{Pi) + N^{P^) + e+{h) + e-{k2). (2.1) 

The invariant eight-fold cross section is 

(fa = , ^ --77^] Y] \T\^ ^dsuds^ (2.2) 
2./A(s, m2,m2) (27r)8 4^^^ ' n! ^ 



X dR2{Pi + P2^q + P12) dR2{q h + Aja) dR2{Pi2 ^ Pi + P2) , 
where the two-body invariant phase space volume R2 is defined as 

dR2{a + b^c + d) = d^P, d'Pd S^^\Pa + Pb-Pc- Pd) 5{Pl - ml) 6{P^ - ml). (2.3) 

The four-momenta of initial (Pi,P2) and final {P[,P2) nucleons are P — {Ep,P) with 
Ep — y/m?~+P^; an analogous notation is used for the lepton momenta ki^2', fn denotes 
the nucleon mass, while the electron mass can be neglected for the present kinematics. 
The invariant mass of two particles is denoted hereafter as s with s = (Pi -|- P2)^; along 
with this notation for the invariant mass of the virtual photon throughout the paper 
we also use the more familiar notation with = M^. The kinematical factor A is 
X{x'^ , y'^ , z"^) — (x^ — {y + z)'^){x'^ ~ {y ~ ^)^); the factor 1/n! accounts for n identical 
particles in the final state. 



B. Leptonic tensor 



The di-electron production process is considered as decay of a virtual photon produced 
in strong and electromagnetic NN interactions from different elementary reactions, e.g., 
bremsstrahlung, Dalitz decay, vector meson decay etc. [20]. For such a process the general 
expression for the invariant amplitude squared reads 

\T\'-W^,p^\ (2.4) 

where the momentum of the virtual photon is denoted as q — {ki + k2); e is the elementary 
electric charge. The purely electromagnetic decay vertex of the virtual photon is deter- 
mined by the leptonic tensor l^^^ — j^j'^ with the current — u{ki,Si) ^^v{k2,S2), 

spins 

where u and v are the corresponding Dirac bispinors for the outgoing electron and positron. 
The leptonic tensor reads explicitly 

IfMu = 4 {kif,k2^ + ki^k2^ - g^u{ki ■ ^2)) (2.5) 

for unpolarized di-electrons. 

The integral over the leptonic phase space is easily calculated due to its covariance 
and the fact that the only "external" variable on which it can depend is the di-electron 
four- momentum q, 

I lAki. k2, q)dR2{q k, + ^2) = y (-g^,, + ^ j . (2.6) 

Obviously, in virtue of gauge invariance of the electromagnetic tensors, q^l^'^ = qj^'^ — 
q^W^y = q^W^i, = 0, only the first term in the r.h.s. of Eq. (2.6) contributes, so that we 
obtain 



A(g,gl2,M^)A(gl2,m^m') ^ +^ 

dM- 6Ms(47r)5 y "^^^""^""^^^ s\2\{s,m\m^) ^J^' ' ^'''^ 

where dil* and d^\2 are defined in the center of mass of initial and final nucleons, respec- 
tively; a^m stands for the electromagnetic fine structure constant. 



5 



C. Lagrangians and parameters 



The covariant hadronic current is evaluated within a meson-nucleon theory based 
on effective interaction Lagrangians which consist on two parts describing the strong and 
electromagnetic interaction. In our approach, the strong interaction among nucleons is 
mediated by four exchange mesons: scalar (a), pseudoscalar-isovector (tt), and neutral 
vector (a;) and vector- isovector (p) mesons [9, 10, 13, 21]. We adopt the nucleon-nucleon- 
meson (NNM) interaction terms 

jONNa = gaNN^(^^), (2.8) 

Cnn. = -^N^5rr{d^^i.))N, (2.9) 
Cnnp = -QNNp (nj^tN^I) - l^iVa^^rTVa^^j;)) , (2.10) 
Cnn. = -^iviv. [ni,N$^^^ - ^Na,,Nd^^l^ , (2.11) 

where and ^?(m) denote the nucleon and meson fields, respectively, and bold face letters 
stand for isovectors. All couplings with off-mass shell particles are dressed by monopole 
form factors Fm — (A|^ — / (A|f — k\j), where k\j is the four- momentum of a virtual 
meson with mass ^m- The effective parameters and their dependence on the initial en- 
ergy are adjusted to experimental data on NN scattering at the considered intermediate 
energies [13, 22]. 

D. Gauge invariance 

The form of the cross section Eq. (2.7) exploits essentially the gauge invariance of 
hadronic and leptonic tensors. This implies that in elaborating models for the reaction 
(2.1) with effective Lagrangians, particular attention must be devoted to the gauge invari- 
ance of the computed currents with the mandatory condition q^J^ — Q. In our approach, 
i.e., in the one-boson exchange approximation (OBE) for the strong A'^A^ interaction and 
one-photon exchange for the electromagnetic production of e+e^, the current is de- 
termined by diagrams of two types: (i) the ones which describe the creation of a virtual 
photon with > as pure nucleon bremsstrahlung as depicted in Fig. 1 and (ii) in case 
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a) b) c) d) 

FIG. 1: Bremsstrahlung diagrams for the process Ni + N2 ^ N[ + N2 + e'^ + e'^ . Fat lines are 
for resonances considered in section III. 

of exchange of charged mesons, the emission of a virtual photon (7*) from internal meson 
lines, see Fig. 2a. For these diagrams the gauge invariance is tightly connected with the 
two-body Ward-Takahashi (WT) identity 

q,T^{p',p) = {^S-\p') - S-\p)) , (2.12) 

where F'^ denotes the electromagnetic vertex and S{p) is the (full) propagator of the 
respective particle. It is straightforward to show that, if (2.12) is to be fulfilled, then pair- 
wise two diagrams with exchange of neutral mesons and pre-emission and post-emission 
of 7* (cf. Figs, la) and b)) cancel each other, hence ensuring g'^J^ = 0, i.e., current con- 
servation. This is also true after dressing the vertices with phenomenological form factors. 
However, in case of charged meson exchange the WT identity is not any more automati- 
cally fulfilled. This is because the nucleon momenta are interchanged and, consequently, 
the "right" and "left" internal nucleon propagators are defined for different momenta of 
the exchanged meson. 

In order to restore the gauge invariance on this level one must consider additional 
diagrams with emission of the virtual photon by the charged meson exchange (Fig. 2a) 
which exactly compensates the non-zero part of the current divergence, and thus gauge 
invariance is restored. This holds true for bar NNM vertices without cut-off form factors. 
Inclusion of additional form factors again leads to non-conserved currents. There are 
several prescriptions of how to preserve gauge invariance within effective theories with 
cut-off form factors. 
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FIG. 2: Contribution of meson exchange currents (a) and seagull terms (b, c) to the process 
-^1 + -^2 N[ + N2 + e'^e~ , where A^i and N2 stand for protons and N2 and denote neutrons. 

The main idea of these prescriptions is to consider the cut-off form factors as phe- 
nomenological part of the self-energy corrections to the corresponding propagators; the 
full propagators are to be treated as the bare ones multiplied at both ends by a form 
factor [23]. Then, the full propagator, e.g. for mesons, can be defined as 



where n(fc) is the self-energy correction (see Fig. 3). 

n(A;) FM{k) Fm{Iz) 



FIG. 3: Graphical illustration of the cut-off form factor as self-energy corrections to the full 
propagator [23]. 

In the simplest case, for mesonic vertices with pseudoscalar couplings, the bare mesonic 
vertex F^^ = {ki^ + ^2^^) receives an additional factor [8, 24, 25] becoming 

^ j^^^ _|_ -J (Aj/ ~ kf) {A\.j — fcl) / ^ kj — ^ kl — \ ^2 

The above prescriptions for restoration of the gauge invariance in pn collisions are 
valid only for pion exchange diagrams with the interaction vertices independent of the 
momentum k of the exchanged meson, i.e., solely for the case of pseudo-scalar irNN 
coupling. The presence of field derivatives in the interaction Lagrangian, e.g. the case 
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for pseudo- vector ttNN coupling or for vector mesons, Eqs. (2.9) and (2.10), requires 
a more refined treatment of the gauge invariance. In the simplest case, besides the WT 
identity condition (with full propagators) for the diagram 2a, the gauge invariance requires 
an introduction of covariant derivatives, i.e. the replacement of the partial derivatives, 
including the NNM vertices, by a covariant form (minimal coupling). Such a procedure 
generates another kind of Feynman diagrams with contact terms, i.e., vertices with four 
lines, known also as KroU-Rudermann [14] or seagull like diagrams, see Figs. 2b and c. 
We include therefore in our calculations these diagrams by the corresponding interaction 
Lagrangian 

jC-NN.-y - -^^iV757M^(T#M)A^ (2.15) 

with electromagnetic four-potential and charge operator e of the pion. Analogously 
for the pNN coupling one has to replace 

2m ^ 2m 2m (''^ ^ ' 

Gauge invariance is henceforth ensured. It should be stressed that as far as the tensor 
part of the pNN Lagrangian (see Eq. (2.10)) is accounted for, the prescription (2.16) 
must be mandatorily apphed, regardless of the choice of ttA^A^ coupling. This implies 
that calculations with pseudo-scalar couplings for the pion-nuclcon vertex violates gauge 
invariance due to p meson exchange. Our numerical calculations show that the effect of 
such p exchange seagull type diagrams vary from 10% at low di-electron invariant masses 
up to 35% at the kinematical limit. 

All electromagnetic A'"A'"7 vertices correspond to the interaction Lagrangian 

- -e + eKN [^^^'') N (2.17) 

with the field strength tensor J^^^, — d^^A^ — d^A,^, and k as the anomalous magnetic 
moment of the nucleon {k, — 1.793 for protons and k — —1.913 for neutrons). 



E. Results for bremsstrahlung 

The OBE parameters and their energy dependence have been taken as in Ref. [8, 
13]. Figure 4 exhibits results of our calculations of the invariant- mass distribution of di- 
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electrons in pp and pn collisions from bremsstrahlung processes in Figs. 1 and 2 (nucleons 
only) at two values of the kinetic energy, 1.04 GeV and 1.25 GeV as relevant for DLS 
[2] and HADES [3] measurements. In our actual calculations we include, besides the 
mentioned four exchange mesons tt, cr, p and co also a "counter term" simulating a heavy 
axial vector-isovector meson, with the goal to cancel singularities of the pion potential at 
the origin [13]. The dotted lines in Fig. 4 depict the cross section in pn collision, while the 
solid lines stand for results of pp reactions. It is seen that thepn cross section is by a factor 
5 — 6 larger than the pp cross section. (The situation for real photon emission is similar: 
The pn channel has a significant contribution, while, due to a destructive interference, 
the pp channel is much weaker and is often neglected [26].) This is due to isospin effects 
for the charged tt and p mesons and different interference effects in pp and pn channels. 
In pn reactions, additional contributions stem from the emission off charged exchange 
mesons and from corresponding seagull type diagrams. Note that in both channels, pp 
and pn, numerical tests of gauge invariance can serve as additional check of the code. 
The meson exchange diagrams together with contact terms amplify the contribution of 
pure nucleonic currents; their contribution is of the order of ~ 40% at low values of the 
invariant di-electron masses and increases up to a factor 2 — 3 at higher invariant masses. 




M[GeV] 

FIG. 4: Contribution of the bremsstrahlung diagrams in Figs. 1 and 2 (without nucleon reso- 
nances) to the e+e~ invariant mass distribution at two kinetic energies (left: 1.04 GeV, right: 
1.25 GeV). Sohd Unes correspond to pp reactions; dashed lines are for pn reactions. 



10 



III. RESONANCES 



Intermediate baryon resonances play an important role in di-electron production in 
NN collisions at beam energies in the 1 - 2 GeV region [13, 16-20, 22, 24, 27]. The main 
contribution to the cross section stems from the A isobar [13]. Also, the low- lying nucleon 
resonances such as A^*(1440), A^*(1520) and A^*(1535) contribute to the cross section. We 
are going to investigate separately each of these resonances represented in Fig. 1 by fat 
lines. 



A. ^33(1232) 

Since the isospin of the A is 3/2 only the isovector mesons tt and p couple to nucleons 
and A's. The form of the effective AN interaction was thoroughly investigated in hter- 
ature in connection with NN scattering [21, 28, 29], pion photo- and electroproduction 

[30-34] . The effective Lagrangians of the ANM interactions read [28-30] 

Can. = ^ {^1 T 9«*(,)7V) + h.c, (3.1) 
jCanp = ^ (^A.T {a/^^f^) - cP^l^ } 7^757V) + h.c. (3.2) 

Ulp 

with fANn = 2.13 GeV and Janp = 7.14 GeV [24]. The ANM vertices are dressed by 
cut-off form factors 



A2 _ i-2 



(3.3) 



where Kan-k = 1-4214 GeV and Aa^v^ = 2.273 GeV [24]. The symbol T in Eqs. (3.1) 
and (3.2) stands for the isospin transition matrix, and ^'a denotes the field describing the 
A. Particles with higher spins (s > 1) are treated usually within the Rarita-Schwinger 
formalism in accordance with which the A propagator has the form 

5f (P, ^a) = llf±^ Pf^p^ ruA) (3.4) 



with the spin projection operator Pa (p) defined as 

2 

\p, niA) = -5"^ + ^7"7^ + + ^ (7 V - 7V) • (3.5) 
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In addition, to take into account the widths of A, the mass in the denominator of the 
scalar part of the propagator is modified as ttia — > ttia — ir^/'^- For the kinematics 
considered here, the "mass" of the intermediate A can be rather far from its pole 
value. The width, as a function of p^, is calculated as a sum of partial widths being 
dominated by the one-pion (A —>■ Ntt) and two-pion (A — >• Np — > Nnn) decay channels 
[35]. 

The general form of the AA'^7 coupling satisfying gauge invariance can be written as 
[32, 36-38] 

^^^^'^Qaf^izsh^TsNd^T"^ + h.c, (3.6) 



4m^ 

Qa^^iz) = ga^,+ [z + ^{l+ 4^)^]7a7M> (3-7) 

where A is a constant reflecting the invariance of the free A Lagrangian with respect 
to point transformations which, according to common practice, is taken A = —1. The 
other parameters, ,21,-22 and Z3 are also connected with point transformations; however, 
they characterize the off-mass shell A resonance and remain unconstrained. The two 
coupling constants gi and g2 in Eq. (3.6) can be estimated from experimental data with 
real photons (g^ = 0) by evaluating the helicity amphtudes [39] of the process A — > 7A'" 
obtained within the Lagrangian Eq. (3.6). One finds explicitly 



egr (T3) = -2x/2^^, IV^A?^^ + A^^''] , (3.8) 
m + rriA V Q ^ 2 J 



e^2(T3) = -4%/2(^^) mm A 



2 mA 2 



(3.9) 



where |q*| is the three-momentum of the photon in the A center-of-mass system, and 
the helicity amplitudes are normalized in such a way that they correspond to 

values hsted by Particle Data Group (PDG) [40]. It should be pointed out that both 
coupling constants 5^1 2 are quite sensitive to the values of the helicity amplitudes, which, 

according to PDG, vary in rather large intervals: — Af-^*^ ~ 0.128 ■■■ 0.145, — Af^"-^ ~ 

2 2 

0.243 ••• 0.261, [40] providing the uncertainties gi ~ 4.5 •••5.5 and g2 ~ 4.5 •••8.5. In 
principle, the computed electromagnetic width of A can serve as an additional constraint 
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for the coupling constants. However it turns out that it is less sensitive to the actual 
choice of gi^2 yielding satisfactorily good results for different sets of gi^2 obtained from 
hehcity amplitudes. Another uncertainty (~ 20%) in Eqs. (3.8) and (3.9) follows from 



the normalization of the isospin transition matrix (T3) which can be chosen either 
or 1 (see discussion appendix A in [8]). The remaining constants in Eq. (3.6), -21,2,3 and 
gs, do not contribute in processes with on-mass shell particles and cannot be directly 
related to data. As a matter of fact, they are considered as free fitting parameters [36]. 
In the present calculation we use gi — 5.478, Zi — 0.05, g2 — 7.611, Z2 — 1.499, 513 = 7.0 
and Z2 — being consistent with the available experimental data on helicity amplitudes, 
electromagnetic decay width and pion photoproduction (cf. Ref. [36]). 
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FIG. 5: Invariant mass distribution of e+e^ pairs in proton-proton (upper row) and proton- 
neutron (lower row) collisions. The dashed (dotted) curves depict the contribution of diagrams 
with bremsstrahlung from jAN {'jNN) vertices. In pn channels, the meson exchange and 
seagull type diagram are accounted for as well. The solid lines are the results of calculations of 
the cross section as coherent sum of all nucleon and A contributions. 



With these parameters we calculate the invariant mass distribution of di-electrons 
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produced in pp and pn collisions. Figure 5 exhibits the mass distribution at two values of 
the kinetic energy, 1.04 and 1.25 GeV. The dotted hues depict the contribution from pure 
bremsstrahlung processes from the nucleon lines including meson exchange and seagull 
diagrams for pn reactions (cf. Fig. 4), while the dashed lines are the contributions of 
the A. The solid lines exhibit the coherent sum of nucleon and A contributions. The 
normalization of the isospin transition matrix has been chosen as 1. It can be seen that 
the A contribution dominates in the whole kinematic range , except for the region near 
the kinematical limit, where the nucleon current contributions become comparable with A 
contributions. It should be noted that, since at considered energies the off-mass shellness 
of A is not large, the contribution of off-mass shell parameters Zi is rather small, except 
in the region near the kinematical limits where the behavior of the cross section is slightly 
modified. However, since the electromagnetic coupling constants Qi have been adjusted to 
experimental data together with off- mass shell parameters [36], we keep Zi as in Ref. [36] 
for the sake of consistency. Also note that the couphng constants gi^2 have been fixed 
at the photon point, i.e., at — 0, while in our case the virtual photon is massive. In 
principle, one may introduce a dependence of the coupling strengths in form of some 
phenomenological form factors which avoid an unphysical behavior at large virtuality of 
A, i.e., at large g^. 

To have an estimate of the role of these parameters for off-mass shell A's, it is instruc- 
tive to investigate the invariant mass distribution in the Dalitz decay of an off-mass shell 
particle with all quantum numbers as the A but with different mass ^ m\. Such a 
quantity is often used in two-step models when calculating di-electrons from NN colli- 
sions [4, 16, 22], where in a first step a A like particle is created, then in the second step it 
decays into a nucleon and a di-electron. There are essentially two options for a treatment 
of such a process: 

(1) Everything is on- mass shell, i.e. in the A Lagrangian, and in the positive energy 
projector operators (p-|-mA), Eq. (3.4), and in the spin projection operator P|(p, ttia), 
Eq. (3.5), one takes the mass parameter as mA =^ rnx- This means that the produced 
on-mass shell particle with a mass mx 7^ is nevertheless described by a A Lagrangian 
with the same coupling constants. This leads merely to a shift in masses in the ex- 
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pression for the invariant mass distribution and to a corresponding enlargement of the 
phase space volume. Evidently, since in this case — m\ and due to the relation 
7q,(p + mx)Pt^{p, TUx) — ip -\- '>nx)Pt^{p, TTix) 7a = 0, the dependence on Zi drops out 

2 2 

in such a treatment (see Eq. (3.7)). 

(2) The A is considered off-mass shell, i.e., in the Lagrangian and projection operators 
(3.5) and (3.4) one keeps the mass parameter as rriA, while p^ ^ m\ (see [16]). In this 
case the positive energy projection operator (p + m/^) does not commute with the spin 
projection operator P| (p, ttia) causing problems in the treatment of the Rarita-Schwinger 
propagator for off-mass shell particles (see discussion in Ref. [8]). Obviously, the off-mass 
shell parameters Zi can now contribute. The mass distribution of the Dalitz decay with 
p^ 7^ m\ but with Ps{p,mA) is called the "off-mass shell" Dalitz decay. 

By using Eq. (2.6), the invariant- mass distribution from the decay of a A-like particle 
into a e+e" pair with invariant mass M can be presented in the following form: 

where |pjv| is the momentum of the outgoing nucleon in the A center of mass system, and 
the amplitude of a virtual photon production Ax — > ^^"7* is defined as 



K^" = TV 



ipN + m) G^'a {p + mx)Pfip,mx)G^,p 

2 



(3.11) 



with 



75- 

(3.12) 

In Fig. 6 (left panel), the mass distribution (3.10) is presented for the case of on-mass 
shell A-hke particles. It can be seen that due to larger phase space volume, the mass 
distribution of heavier particles is much larger than the distribution for real on-mass 
shell A. 

In Fig. 6 (right panel), a comparison of results of calculations of the off mass shell 
distribution arc presented for mx = 1-85 GcV. The solid line is for the "on- mass shell" 
case (1), while the off- mass shell results for case (2) are presented by dashed lines. It can be 
seen that the off-mass shell calculations considerably differ from the commonly accepted 
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M [GeV] M [GeV] 

FIG. 6: Invariant mass distribution according to Eq. (3.10) for the decay of a A-like particle 
mx — Ne^e~. Left panel: the mass parameters mx = 1.232 GeV (lower curves) and at 
mx = 1.850 GeV (upper curves) are taken on-mass shell, i.e., case (1) with = m\. The solid 
and dashed lines correspond to isospin normalization (T3) to 1 and to respectively. Right 

panel: The on-mass shell mass distribution (solid line) versus off-mass shell calculations in case 

(2), i.e., with the spin projection operator Ps (p, mx), Eq. (3.5), and the electromagnetic vertices 

2 

Eq. (3.12) calculated at 7^ m\ (dashed line). In all calculations the coupling constants in 
Eq. (3.12) are g\ = 5.478, g2 = 7.611 and 53 = 7.0. For the right panel the off-mass shell 
parameters have been taken as Zi = —0.5. 

results with A on-mass shell. This demonstrates that in two-step models, besides the 
traditional approximations, there are additional uncertainties in the treatment the Dalitz 
decay of the off-mass shell A at large values of invariant mass. It is clear that, if in 
two-step models off-mass shell calculations are employed, additional form factors must be 
considered to suppress such an increase of the mass distribution at large virtuality of the 
decaying A. 

B. Spin-i resonances Pii(1440) and 511(1535) 

Contrarily to the A isobars, the isospin- 1 nucleon resonances can couple not only 
with isospin-1 mesons, but also a, r] and lu mesons may contribute. The corresponding 
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Lagragians are 



'(±) 

■'NN*ps 



dNN' 



ps 



(±) 

NN*V 



rriN* ± m 

gNN*v 
2{mN* + rn) 



N* 



75 
1 



1 

75 



a^,V^"'{x)N + h.c., 



1 

75 



A^%) + /i.e. 



(3.13) 
(3.14) 
(3.15) 



with the abbreviations ps = tt or rj, ^{ps) = {T^{n)) or $(^)(a;), V = V((^) or V{Tp), and 
_ Ql3ya _ Qayp ^ (The additionally introduced r] exchange employs for the NNr] 
interaction C^n,, = -^^N-f5-f^'{d^^^r,))N with f NNr, = 1-79 from Ref. [21].) At a first 



glance, the consideration of Lagrangians (3.13)-(3.15) involves into the calculations ad- 
ditional free parameters. However, the effective constants are not completely free and 
indeed they can be estimated from independent experimental data. The coupling con- 
stants of the resonances with the pseudo-scalar meson can be found from the calculations 
of the experimentally known [40] partial width of the corresponding decay 

X |p' 



J- N'^Nps 



2 

Qnn* 



(3.16) 



where T = 3 for pions and X — 1 otherwise. This yields glf^'n — 6.54, g'^jv*?? — O-^ the 
Roper resonance A'"(1440) and 9%^*^^ — 1.25, glfN*r) — 2-02 for the A'^(1535) resonance. 

Estimates of coupling constants with a mesons, for which direct experimental data 
are not available, are more involved. One of the method is to calculate the decay of the 
considered resonances into two pions in a specific final state, N{2TT)glPQ. For low-lying 
resonances such a decay can be treated as a process with intermediate excitation of the 
a meson with its subsequent decay into two pions [41], as depicted in Fig. 7. 

A direct calculation of the diagram Fig. 7 results in 



JV*^iV(27r)^=o 



9nn* (t 
47r 



nif^* —m 



|p*|m^ [{rriN* ± m)^ - V 



2mTv 



\e - (m, - in 



212 



^2 _ 4^2 



^,(3.17) 



where the dependence of the width cr — 27r on the mass ^ of the intermediate meson is 
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N* 



FIG. 7: Diagram for the decay of a resonance N* in to two pions in the state = 0, / = with 
an intermediate a meson. 



computed from the Lagrangian 



-^^aTTTT = S'tTTTTT^ (*(7r)*{7r)) "^"((7) (3.18) 



2 

with the coupUng constant (^o-tttt found from the total decay width of the a meson into two 
pions. In the present calculation we adopt m„ = 500 MeV and ro-^27r = 250 MeV being 
consistent with the recent analysis [42]. With these parameters we obtain gNN*cr = 2.1 
for the Roper resonance and gNN*a = 3.8 for A^(1535) (see also [43]). The propagator of 
the off-mass shell resonance, 7^ fnj^,, is augmented by a form factor of the form 



with A = 1.2 GeV in both cases, for Roper and A^(1535) resonances. The finite widths 
of the intermediate resonances are taken into account as usually: adding to the mass 
parameter in the propagator an imaginary part as miv* — m^v* — ir*°*/2- In our case, 
the resonances are off mass shell and their total widths depend on the invariant masses 
of the resonance. Such a dependence can be parameterized as [36] 

T'°\mx) = E Ti{mx)Fx{mx), (3.20) 

i 

where the sum runs over all possible partial decay channels, and the cut-off form factors 
Fx{fnx) suppress an unphysical increase of the width with increasing mx- 

Fximx) = , pL,,x» , (3.21) 



,P*(mjv*). 

where P* is the nucleon momentum in the resonance center of mass system, a = 3 and 
a = 2 for the Roper and A^(1535) resonances, respectively. For the Roper resonance 
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there are two main decay channels, A^(1440) Ntt (branching ratio ~ 60 — 70%) and 
A^(1440) Att (branching ratio ~ 20 - 30%) [40], while A^(1535) decays mainly either 
into a nucleon and a pion (branching ratio ~ 55%) or into a nucleon and rj (branching 
ratio ~ 45%). The energy dependence of the partial widths Ti{mx) for two-body decay 
have been calculated using the same Lagrangians (3.13)-(3.15). For the decay of the Roper 
resonance A^* — > Ntttt via A we employ an effective Lagrangian of the form 

Cn^a. = ^^^*-T9,*(,)7V* + h.c. (3.22) 

Note that by calculating the energy dependence of the widths from the Lagrangian (3.22) 
a knowledge of the coupling constant Jn^Att is not necessary. Also note that, in spite of 
the mass of the Roper resonance which is only slightly above the kinematical limit, the 
probability to decay into A and tt is relatively large. This is due to the large total width 
of the A resonance which correspondingly spreads the mass around the pole position, 
enlarging therefore the phase space of the decay channel. In calculating the partial width 
N* ttA we adopt a Gaussian distribution of the mass of the A resonance /(ttia) = 
((TA'\/27r)~^ exp[— (rrtA — "^A)^/(2cr|^)] with cta = Ta- The remaining coupling constants 
for the vector mesons have been taken from Ref. [43] (see Tab. I) . 

TABLE I: Coupling constants gNN*M and cut-off parameters A for the effective Lagrangians 
(3.13)-(3.15), computed either directly from the decay widths or taken from Ref. [43]. 



Meson 


Pn(1440) 


Sii(1535) 


A [GeV] 


TT 


6.54 


1.25 


1.2 


V 


0.5 


2.02 


1.2 


a 


2.1 


3.8 


1.2 


P 


-0.57 


- 0.65 


1.2 




- 0.37 


-0.72 


1.2 



The effective Lagrangian for the electromagnetic decay of the resonance into a photon 
and a nucleon has been taken as 

1 



'-'NN*'Y 



eK 



2mR 



N* 



75 



} a^'NF., + h.c, 



(3.23) 
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where, in contrast to the A case, the couphng k is different for proton and neutron vertices 
and can be found from the hehcity amphtudes 

Ay2^-eK\U^, (3.24) 

where for the Roper resonance one has A^^^ — —0.065 GeV~^/^ and A^^^ — 0-04 GeV^^/^, 
while for the A^(1535) we employ A^^^ = 0.09 GeV-V^ and = -0.046 GeV-V2, 
correspondingly. 



C. Spin- 1 nucleon resonances 



The next considered resonance is Di3(1520) with negative parity and spin |. The 



2' 



effective Lagrangian for such resonances is chosen in the same form as for the A with the 
exception that, since the isospin is |, besides the isovector mesons tt and p, the isoscalar 
77, cu and a also can contribute. The effective Lagrangians are as follows 



'-'NN*ps — ^N*\-^) ' 

Tflpg 



75 



C 



(±) 

NN*V 



'(±) 
■'NN*a 



— I 




{da^(ps))N + h.c., 



> {da^(a))N + h.C.. 



(3.25) 
(3.26) 
(3.27) 



Note that in choosing the relative phase for the a meson an imaginary unit i must be ex- 
plicitly displayed. In principle, to synchronize the relative phases of different Lagrangians 
one may compute the corresponding amplitude in a fully coplanar kinematics. Then such 
an amplitude, in tree level calculations, must be either purely real or purely imaginary. 
The coupling constants and the cut-off parameter A = 0.8 GeV have been taken from 
Ref. [13]. The propagator is chosen in the form (3.4) with the resonance mass augmented 
by the total decay width. The total width of the resonance A'"(1520) is calculated by 
Eq. (3.20) for which three decay channels have been taken into account, A^(1520) Ntt 
(branching ratio ~ 50%), Ar(1520) Np (branching ratio ~ 25%) and iV(1520) NA 
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(branching ratio ~ 25%). For the decay in to Np and A'"A the partial widths have been 
calculated again by adopting Gaussian distributions of the mass of p and A around their 
pole values. 

The electromagnetic part of the Lagrangian has the same form as for the A, see 
Eq. (3.6), except for the isospin transition matrix T, i.e. 



As in case of A, the two coupling constants g^i 2 can be obtained from the helicity am- 
plitudes, separately for proton and neutron. The remaining constants g^, -21,-22 and 2:3, 
being completely free parameters, are to be found by fitting experimental data. There 
exist several sets in the literature equally well describing the corresponding data [36]. 
In our calculations we have chosen gi — 3.004, g2 — 3.047, g^ — Q for the proton and 
gi — —0.068, g2 — 1.265, g's = for the neutron. The off-mass shell parameters zi — 1.39 
and Z2 = 0.267 have been also taken from [36] (since fifs = 0, the off mass shell parameter 
^3 is irrelevant here). 

As in case of A isobar we calculate the mass distribution (3.10) of the Dalitz decay 
of a particle with quantum numbers of A'^(1520) at different masses. This distribution 
is also frequently used in calculations by two-step models [4]. Figure 8 illustrates the 
behavior of the invariant mass distribution in the Dalitz decay of a Dia-like resonance 
at different masses (left and middle panel) for the calculations of on-mass and off-mass 
shell decays. It can be seen that, as in the case of A isobar (cf. Fig. 6), the heavier 
masses result in broader distributions. The virtuality of the decaying resonance leads to 
an increased broadening distribution. Also, since the neutron couplings 511^2 are much 
smaller than the proton ones, the effect of virtuality is less pronounced in the neutron 
case. In the right panel of Fig. 8 we present a direct comparison of the contribution of the 
A^(1520) resonance for proton and neutron decays. At larger di-electron invariant masses 
the relative contribution of the A^(1520) in pn collisions is smaller than in pp processes. 




(3.28) 
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FIG. 8: Invariant mass distribution according to Eq. (3.10) for the decay of a Ar(1520)-like 
particle mx — > Ne^e~ into a proton (left panel) and a neutron (middle panel) at two values of 
the mass parameter mx = 1.520 GeV and mx = 1.850 GeV. For the case of mx = 1.850 GeV 
two different definitions of dV^^^'^^^ /dM have been used corresponding to cases (1) and (2) 
in subsection III A (see text there). The solid and dotted lines correspond to the case (1), where 
the resonance is on mass shell; the dot-dashed lines are results of calculations for case (2) with 
the resonance off mass shell, i.e., in the spin projection operator P3{p,mx), Eq. (3.5), and the 
electromagnetic vertices (3.12) the mass parameter is mx = 1.520 GeV, while = 1.850 GeV. 
In the right panel, the comparison of the mass distribution for proton and neutron vertices 
reflects the relative contribution of the iV(1520) resonance in pp and pn di-electron production 
(cf. solid curves in the left and middle panels). 

With the above parameters we calculate the contribution of the mentioned baryon 
resonances in the invariant mass distribution of di-clcctrons from pp and pn collisions. 
In Fig. 9 we present the comparison of the nuclcon bremsstrahlung contribution with 
individual contributions from each of the considered resonances in pp (upper row) and pn 
(lower row) collisions at two kinetic energies. 

The role of the Roper resonance (dash-dot-dot lines) is neghgibly small in all cases. 
Also a small contribution comes from A'"(1535) (dash-dotted hnes) which becomes of the 
same order as the nucleon bremsstrahlung only at the kinematical limit. A more significant 
contribution stems from the A^(1520) resonance, which becomes competitive with nucleon 
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M [GeV] M [GeV] 

FIG. 9: Invariant-mass distribution of e'^e~ in pp (top) and pn (bottom) collisions. The dotted 
curves depict the contribution of diagrams with bremsstrahlung from ^NN vertices, cf. Fig. 4. 
Dash-dot-dot lines: contribution of Roper-resonance, dash-dotted lines: A'^(1535), dashed lines: 
iV(1520). Solid curves are for all resonances including the contribution from A. 

bremsstrahlung already at di-electron invariant mass M > 0.45 GeV. The sohd hne is the 
coherent sum of all the resonances, including A isobars. In these calculations and in what 
follows, the isospin transition matrix has been normalized to ^2/3. The contribution 
of A'"(1535) becomes more pronounced at di-electron invariant masses corresponding to 
the pole position of the resonance mass. Note that in our calculations the initial state 
interaction has been taken into account by imposing a energy dependence of the effective 
parameters as suggested in Ref. [13]. In principle, one could account for the effects of 
initial state interaction explicitly, as proposed in [44]. An analysis performed in [15] 
shows that at the considered energies the two methods of accounting for the initial state 
effects provide similar results. The effects of final state interaction (FSI) in the considered 
reactions have been investigated in Ref. [8] . It has been found that FSI corrections depend 
on the relative momentum of the outgoing nucleons, becoming significant at low momenta. 
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i.e. at the kinematical limit of the di-electron invariant mass. At low and intermediate 
values of the di-electron mass FSI effects are small. 




M [GeV] M [GeV] 

FIG. 10: Invariant mass distribution of e'^e~ pairs in pp and pn collisions at two kinetic ener- 
gies (left: 1.04 GeV, right: 1.25 GeV). The dotted (dash-dotted) hnes depict contributions of 
bremsstrahlung diagrams without resonances in pp (pn) collisions. The solid lines are the re- 
sults of the coherent sum of all the diagrams, including bremsstrahlung and contributions from 
P33(1232), Pii(1440), £)i3(1520) and 5ii(1535) resonances. 

Eventually, the total cross section with accounting for all resonances and nucleon 
bremsstrahlung is presented in Fig. 10 by solid lines. The dotted and dash-dotted curves 
exhibit the contribution from the nucleon bremsstrahlung solely in pp and pn collisions 
respectively. One concludes from Fig. 10 that the contribution of baryon resonances 
dominates the cross section at the considered energies. 

It is worth emphasizing that in the present calculations the quantum mechanics in- 
terference effects play an important role in the total cross section, essentially reducing 
the cross section in comparison to a incoherent sum of different contributions. This is 
illustrated in Fig. 11, where results of a coherent summation of Feynman diagrams are pre- 
sented and compared with the incoherent sum of separate contributions of bremsstrahlung 
and P33(1232), Pii(1440), Dis{1520) and ,Sii(1535) resonances. It is seen that in both 
cases, pn and pp collisions, the interference effects become significant at higher values of 
the di-electron invariant mass and reduce the cross section by a factor of about 2 — 2.5. 
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FIG. 11: Invariant mass distribution of e~^e~ in pp and pn collisions as a coherent sum of 
the considered Feynman diagrams (solid lines) vs. an incoherent summation (dashed lines) of 
separate contributions from bremsstrahlung, P33(1232), Pii(1440), Di3(1520) and 511(1525), 
respectively. 

Experimentally, information on the di-electron production from pn collisions may be 
extracted from the tagged neutrons in Dp Psp^p e'^e" reactions by exploiting the so- 
called spectator mechanism. As discussed in Ref . [8] , if the spectator proton psp is detected 
in the very forward direction with about half of the momentum of the incident deuteron, 
then with a high probability the reaction occurred at the neutron, and the proton from 
the deuteron remains as a spectator. In such a case, one may extract the pn sub-reaction 
at the same energy as the detected proton. To reduce the experimental errors one may 
measure the ratio (Jpn/'^pp in such experiments. However, even such a ratio may remain 
rather sensitive to the extraction procedure, namely to the accuracy of determining the 
effective momentum of the tagged active neutron. 

In Fig. 12 we present the ratio o^nl^pp calculated at few different kinetic energies in the 
pn collisions while keeping fixed the kinetic beam energy of 1.25 GeV for the pp reaction. 
Such a ratio emulates roughly the possible Fermi motion effects in the Dp — > Ps^np e^e~ 
subreaction. It can be seen that, for di-electron invariant masses M > 300 MeV, the 
presented ratio is quite sensitive to the effective momenta of the neutron. 
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FIG. 12: Ratio of the e+e" invariant mass distribution in pn and pp collisions at kinetic beam 
energy of 1.25 GeV for the pp reaction and four different energies for the pn reaction: 1.20 GeV 
(dash-dotted line), 1.25 GeV (solid line), 1.30 GeV (dotted line) and 1.35 GeV (dashed line). 

IV. SUMMARY 

In summary we have analyzed various aspects of the di-electron production from the 
bremsstrahlung mechanism and resonance excitations at intermediate energies for the 
exclusive reactions NN — > NN e'^e^ , i.e., for pp ppe^e~ and np — > npe^e~ . To 
calculate the corresponding cross sections we employ an effective meson-nucleon theory 
with parameters adjusted to elastic A^'A'" and inelastic NN NNtt reaction data with 
low-mass baryon resonances included. 

The performed evaluations of bremsstrahlung diagrams can be considered as an esti- 
mate of the background contribution, a detailed knowledge of which is a necessary prereq- 
uisite for understanding di-clcctron production in heavy-ion collisions. Our approach is 
based on covariant evaluations of the corresponding tree level Fcynman diagrams with im- 
plementing phenomenological form factors, with particular attention paid on preserving 
the gauge invariance. It is stressed that, regardless of the choice of the pion-nucleon- 
nucleon coupling, the consideration of seagull type diagrams is inevitable if meson field 
derivatives enter the interaction Lagrangians, say for p mesons. The covariance of the 
approach is ensured by direct calculations of Feynman diagrams. 



1.20 GeV 
•1.25 GeV 
1.30 GeV 
1.35 GeV 
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In accordance with previous results [8, 13] our calculations demonstrate that in the 
region of invariant masses sufficiently far from the vector meson production threshold 
the main contribution to the cross section, in both reactions pp and pn, comes from 
virtual excitations of nucleon resonances. The contribution from the Roper resonance 
is negligibly small in all the considered reactions. The role of the 5*11(1535) resonance 
is also marginal in the whole kinematical region except for values of the invariant mass 
at the kinematical limit. The main contribution to the cross sections comes from spin- 
I resonances, A and A^(1520), the role of the latter increasing with increasing initial 
NN energy. Due to isospin effects and meson exchange diagrams the cross section for 
the reaction pn ^ pn e'^e" is larger than the cross section for pp pp e'^e" by a factor 
1.5—3. Note that because of (i) contributions of the isoscalar a and uj exchange mesons, (ii) 
differences in the electromagnetic coupling in 7p and 'jn systems, (iii) interference effects, 
and (iv) contribution of resonances, the isospin enhancement is not ~ 9, as one could 
naively expect from isospin symmetry considerations. In both reactions, pn — > pn e^e~ 
and pp — > pp e'^e", the bremsstrahlung cross sections exhibit a smooth behavior as a 
function of the di-electron mass. Hence, the bremsstrahlung cross section can indeed be 
considered as background contribution. 

The previous "DLS puzzle", experimentally resolved in [3], seems to be shifted now to 
a "theory puzzle": the preliminary data for the invariant mass spectrum in the reaction 
np — > np e'*"e~, extracted from the tagged subreaction in Dp — > PspUp e'^e~, point to a 
shoulder at intermediate values of the di-electron invariant mass [45] . Such a structure is 
hardly described within the present approach. (In contrast, the use of the phenomenologi- 
cal one-boson exchange model for the exclusive reactions A'"A'" NNM with M = uj, (p,ri 
and 1]' [9-11, 46] successfully describes data and has some prediction power [47]). Thus, 
understanding the elementary channels remains challenging. Finally, it should be em- 
phasized that we consider here the exclusive reaction NN NN e^e~. The inclusive 
channels NN — > NN X e^e~ may be significantly different (cf. [48]). 
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